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We develop a perturbative method of computing spectral singularities of a Schrodinger operator 
defined by a general complex potential that vanishes outside a closed interval. These can be realized 
as zero-width resonances in optical gain media and correspond to a lasing effect that occurs at 
the threshold gain. Their time-reversed copies yield coherent perfect absorption of light that is 
also known as antilasering. We use our general results to establish the exactness of the 71-th order 
perturbation theory for an arbitrary complex potential consisting of n delta-functions, obtain an 
exact expression for the transfer matrix of these potentials, and examine spectral singularities of 
complex barrier potentials of arbitrary shape. In the context of optical spectral singularities, these 
correspond to inhomogeneous gain media. 
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INTRODUCTION 



A spectral singularity is a well-known mathematical concept [l[ with an interesting physical counterpart; it cor- 
responds to a zero- width resonance @- As shown in Refs. [U, Hj], such resonances can be realized in optical systems 
consisting of a gain medium. They give rise to a particular lasing effect that occurs at the threshold gain Q ■ Since the 
publication of [2J there has been a growing interest in the study of the physical applications of spectral singularities 
In particular, it turns out that a time-reversed copy of an optical spectral singularity (OSS) corresponds to 
a coherent perfect absorption of light Q . This is the basic physical phenomenon occurring in an antilaser [ij . 

The investigation of OSSs that is conducted in relies on the assumption that the gain (attenuation) coefficient 
is constant throughout the gain (loss) region(s). This is a simplifying condition that is almost impossible to fulfill in 
practice. In order to have more realistic models displaying an OSS we need to develop computational techniques for 
the cases that the optically active medium is inhomogeneous. The first step in this direction is taken in Q, where a 
semiclassical method of calculating spectral singularities is employed. The purpose of the present paper is to develop 
an alternative method of exploring spectral singularities that is based on perturbation theory. This is especially 
useful, because for most optically active media, the imaginary part of the effective potential, that is responsible for 
the emergence of an OSS, is several orders of magnitude smaller than its real part |10| . 

Consider the time-independent Schodingcr equation 



i)"{x) + v{x)ip{x) = k 2 ip(x), 



(1) 



where ieR and v is a complex- valued potential that vanishes outside the interval [0, 1]. We can express the general 
solution of fl]) as 



A_e lkx + B-e~ lkx for x < 0, 
ip(x) = { A (j)i{x]k) + B (t) 2 (x; k) for < x < 1, 

for x > 1, 



A + e ikx 



B+e -ikx 



(2) 



where A±, B±, Aq, Bq are complex coefficients and <pi(-;k) and 4>2(',k) are solutions of |T|) on the interval [0,1] that 
are determined by the initial conditions 0: 



0i(O;fc) = l, cf>[(0;k) = -ik, <f> 2 (0;k) = l, <f>' 2 (0;k) = 0. 



(3) 



Spectral singularities are given by the real zeros of the M22 entry of the transfer matrix M of the system Q ■ This 



is the 2x2 matrix M that satisfies 



A, 
B 4 



= M 



A- 
B_ 



Demanding that ip is a continuously differentiable function 



(throughout R) and using Q, we find the following expression for the transfer matrix. 

-e- lfc [r 1+ (fc) - 2T 2 +(fc)] e-*T 1+ (fc) 
e lk [IV (k) - 2IV (k)} -e ifc ri_ (fc) 



M = 



2ik 



(4) 
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where 

T j± (k) :=4f j (l]k)±ik<P i (l;k). (5) 

In view of (U]) spectral singularities are the real zeros of the Jost function Tx-i Q- 

In order to develop a perturbative method of computing spectral singularities we consider potentials of the form 

{vq(x) + evi(x) for < x < 1, 
(6) 
otherwise, 

where vq is an exactly solvable potential, e is a real perturbation parameter, and v\ is an arbitrary potential. Our 
aim is to use perturbation theory to compute the solution (f>i(-; k) of the Schrddinger equation ([1]), the Jost function 
ri_, and its real zeros. 

The plan of the paper is as follows. In Section 2 we derive a perturbative series expansion for </>i(-; k) and discuss 
how it can be used to locate the spectral singularities. In Section 3 we use our perturbative method to examine 
potentials involving one or more delta functions. Here we establish the exactness of the perturbation theory and give 
an explicit formula for the transfer matrix of the system. In Section 4, we apply our method to an arbitrary piecewise 
continuous complex barrier potential. In Section 5, we use the results of Section 4 to study the OSS of an infinite 
planar slab gain medium whose gain/loss coefficient varies along the normal direction. Finally, in Section 6 we give a 
summary of our findings and concluding remarks. 



II. PERTURBATIVE CALCULATION OF SPECTRAL SINGULARITIES 



We start our analysis by introducing the differential operator 

d 2 

L:= — -v Q (x) + k 2 , (7) 

and note that </>i(-; k) and </>2(s k) are solutions of the differential equation 

Ltf>(x) = evi(a#(z), (8) 

on the interval [0, 1] that are uniquely determined by the initial conditions ([3]). To construct a perturbative solution 
of ([8|) we insert the ansatz 

oo 

#aO=£>«>(a;)e* (9) 

3=0 

in ((5J and demand that it holds term by term in powers of e. This yields 

L^°\x) = 0, (10) 
L4> {t) {x) = vxixW-^ix), (11) 

where £=1,2,3,---. We can express the general solution of the latter equation in the form 

r ) (x)=a e ^\x)+b e ^\x)+ f G(x,y)v 1 (yW- 1) (y)dy, (12) 

Jo 

where ae,be are constant coefficients, tp^ and 02 are linearly-independent solutions of (TTU)) , and G is the Green's 
function for the operator L that can be expressed as [ll[ : 

Gfeg)= ^wfw-4°'wS°'w. (13) 

4> l i'(v)*f 6/)-4 0) &¥i 0) (v) 

A convenient choice, that determines <j)f^ and <j)^ in a unique manner, is to demand that they also satisfy the 
initial conditions i.e., 

0i O) (O) = l, 4° } '(0) - -ifc, <4 0) (0) = 1, 4 0) '(0)-0. (14) 
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In this case the denominator of the Green's function that coincides with the Wronskian of <p^ and <j>^ takes the 
value ik, and we find 



G(x, y) = ik- 1 4 0) {x)4>T (y) - <$ ] (x)^ (y) 



(0), 



(,(°)/ 



(15) 



Another advantage of imposing (|14[) is that it identifies the zeroth-order term in the perturbative expansion of </>j (-; k) 



with <pf (•; k). Moreover, if we denote the higher order terms in this expansion by <j) ; (-; k) t l (with i > 1), so that 



(0, 



Mx; k) = J2 4>f{x\ k) /, ^ 0) (x; k) = <j>f (x), 



(16) 



then ([3]) and (| 14(1 imply that (^(O; fc) = (f>^ ' (0; fc) = for £ > 1. A direct consequence of this relation is that the 
coefficients ag and bg of (fT2")l vanish, and we can use (fT2"]l to derive the formula: 



<t>f{x t ,k) 



■in • x,„-i)ui(a; m _i), 



(17) 



that holds for all £ > 1 and a;^ G [0, 1]. 

We can use (|T6| to obtain a perturbative expansion for the Jost functions ([5]): 



r J± (fc) = £ r §( fc ) e ' 



t=0 



In view of ([5]), (JTHJ) , and (|T7)) . we have 



vfl{k) = (l)±ik^'(l), 



- Ao)' 



(0), 



rg(*) 



dxt 



dxe 



dxi cj)f\xi) x 



t — 1 

| [G'(l,a;^) ± ifcG(l,X£)] wi(xf) G(x m+ i,x m )wi(a;. m )|, 



m— 1 



(18) 



(19) 



(20) 



where 



G'(x, y) := 8 x G(x, y) = zfc" 1 0<°> {x)4>T (») - 4^' (*X J M 



/.(°) 



(21) 



Next, we recall that in general vq and f i involve a set of complex coupling constants (31, 32, • ' ' > 3n)- These together 
with the perturbation parameter e are the parameters that enter in the expression (1181) for the Jost solutions. Spectral 
singularities are given by the real values of k for which 



ri_(fc) = 0. 



(22) 



This is a complex equation involving n complex variables, 31,32, ■ ■ ■ ,3m and two real variables, e, k. 

If the unperturbed potential vq is real- valued or more generally docs not support a spectral singularity, the emergence 
of spectral singularities is a consequence of the perturbation e v\ . Otherwise the presence of the latter leads to small 
changes in the location of the spectral singularities of vq that are given by the real solutions of 



rfl(k) = o. 



(23) 



In this case, we can again adopt perturbative theory to construct solutions of (|2"2"j) using those of (|2"3"]) . fl2| . Before, 
we explore the details of this construction, we consider a situation where perturbation theory yields the exact solution 
of the problem. 



III. ARRAY OF COMPLEX DELTA-FUNCTION POTENTIALS 

Consider the case that vq (x) = and the perturbation involves an array of Dirac delta- functions EMI: 

n 

vi{x)= y ^i i 5{x-a i ), (24) 

where 31,32, •' ' >$n are complex coupling constants, and a±, a^-, ■ ■ ■ ,a n & re arbitrary positive numbers satisfying 

< ai < a 2 < ■ ■ ■ < a n < 1. (25) 
Because v = 0, Eqs. ([7]l, ([TO ]) . (fT3 ]l , ([X5 |l . and give 

(j>f\x) = e~ lkx , $\x) = cos(fcx), (26) 

G(x,y) = Sm[Hx k ~ y)] , G'(x,y) = cos[k(x-y)}. (27) 

If we substitute (f2"4"]l . (f2l))) . and (|27|l in ([T7| . and use the properties of the delta- function to perform the relevant 
integrals, we find 



cf>f'(x; k)=Y, W MH* ~ ai)}0(x - a,), (28) 

i=l 

where £ > 1, for all ii = 1, 2, • ■ ■ , n and j = 1, 2, 



Zfu ■= k~ e h e J2 ^K) II 3^ sin[fc(a Jm+1 - a im )]^(a Wl - a im ), (29) 

zi,22,- -- ,i^_i=l m= 1 



and 6 stands for the Heaviside step function, 



In view of (f2"5"j) and (|3T))) . the product on the right-hand side of (f2"9"|) vanishes identically, if i m +i < i TO . Hence 

ie l—l 

Zf^k) = 6{n - 1) k-% £ ^ 0) (°*i) IT 3<™ sin [ fc («^ + i - <01. (31) 

u<*2<-"<^— 1=1 m=l 

and as a result <pj(x\ k) = for all £ > n. This proves the following theorem. 

Theorem: For a point interaction consisting of n delta-functions with arbitrary centers and possibly complex 
coupling constants, the n-th order perturbation theory is exact. 

Next, we compute the Jost functions (US]). Using ([2^|). (gg) , and (gTJ) in (JTSJ) and (|20]). we find 

±1 

ti«2<— <»< = ! P=l m=l 



r™(i) = 0, r^(fc) = -2i e - ifc , r^(fc) = ±^e ±ifc , (32) 

/ 1 1 \ i-l n i £-1 

rg(fc) = 6{n-£) 53 II 3^ IT [l-e^ fc K„ +1 -^)j , (33 ) 



where £ > 1, and for all i\ = 1, 2, ■ 



n <li+ -e^ 1 " 2 ^), e- <fc , Oi l2 ± := ^ (l + e™^) e ±ik . (34) 
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Substituting (f26|) in this relation, using the result in Q, and setting e = 1, we obtain the following expressions for 
the entries of the transfer matrix corresponding to the potential 



t e-i 



Mia = 



M 21 



M 22 = 



-2ik(a-i i 



l-e 



-2ik(ai 



e runt 1 < 

£=1 i 1 <i 2 <-"<i^=l p=l m— 1 

n n £ £-1 

£=1 ii<i2<---<if = l p=l m=l 

n n ^ £-1 

Ec-^r* e ^ n *. n I 1 - e2ife(a>m+i ai 

£—1 i 1 <i 2 <...<j £ =l p=l m— 1 

i+E(-^)-^ e runt 1 



2ik(a,i , . - 
e v m + 1 



(35) 
(36) 
(37) 
(38) 



ii <»2<---<if = 1 p=i m=l 



It is instructive to examine the cases n — 1,2. 

For n = 1, we have vi = 3i<5(x — ai), and (f35|) - (|38|) give 



*3i e 



-2ika\ 



2k 



M 21 = 



«3i e 



2ika\ 



2k 



m 22 = l + Mi, 
2fc 



(39) 



This agrees with the results of II 711 . In particular, a sp ectral singularity occurs for imaginary values of 31 and is given 
by k = — iji/2, as envisaged in |l5| and shown in jig . 

For n = 2, we have v\ = i\8{x — a\) + ^S(x — a 2 ), and (|35|) - (f38|) give 



Mu 
M 12 
M 2i 
M 22 



^(3i +32) 3i32[l-e- 2 * fc ^- ai )] 
2fc 4fc 2 



2k 



4fc 2 



■^ ie 2ifc Ql +32e 2ifca 2 ) ^^(^fca! _ e 2ifea 2 ) 
2fc : 

i(3i+3 2 ) 3i32 

2k 4fc 2 



(40) 
(41) 
(42) 
(43) 



Spectral singularities are therefore given by the real values of k for which the right-hand side of (|43p vanishes. Again 
this is in complete agreement with the results of [l6| . 

For the cases that := i/(n+ 1) and 31 = g 2 = • • • = 3„, the potential « is locally periodic, and we can use the 
results of [14[ to compute the transfer matrix of the system. We have checked by explicit calculation for small values 
of 71 that (|3"5|) - (f3"5)) give the same expression for the transfer matrix as the one derived in [l4[ ■ 



IV. COMPLEX BARRIER POTENTIALS 

Consider the potentials of the form © where for all x € [0, 1], 

va{x):=ii, vi(x) :=i 2 f{x), (44) 

31 and 3 2 are complex coupling constants, / : [0, 1] — > C is an integrable function satisfying J 1 dx\f(x)\ < 1, and 
|e3 2 | <C 1 3 1 1 . As explained in Ref. Q, these potentials appear in the study of the OSS of an infinite planar slab gain 
medium with gain coefficient changing along the normal direction to the slab. 

In order to determine the spectral singularities of the potentials given by © and (|4"4"|) , we first use ([7]), (fTDj) . (|14l) . 
(fT5)) . and (|2"Tj) to compute 

(^^(x) = cos(nfca;) — in -1 sin(nfca;), <j)^\x) = cos(nfca;), (45) 
G(x,y) = (nk)- 1 sm[nk(x-y)}, G'{x,y) - ikG(x,y) = 0<°> (x - y), (46) 
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l-§- (47) 



where 



Clearly, n = marks a singularity of our construction that we will avoid. If we substitute ([45)1 in ([19)) . we find 

ri°2(fc) = F (n,k) := -n^k [(n 2 + 1) sin(nfe) + 2mcos(nfc)] 

(48) 



k(n+l) 2 e 
2in 



2 'ink 



x 2 

n — 1 



n+1 



This is consistent with the results of Ref. because it implies that for e = 0, that corresponds to a constant complex 
barrier potential, the spectral singularities are determined by the equation: 

e -2ink_f^_zl\ =Q (4g) 

Similarly, using (|20]). ([4"5j) and (gg]), we obtain for all i > 1, 

T ( (l{k)=iiF l {n,k) (50) 

where 

"1 /-cc/ pxo £—1 



F e (n,k) :-- 



(nA;) 1 "^ / dx t / d&t-\--- \ dx 1 ^(n,k,x 1 ) TJ sm[nk(x m+ i - x m )]T\ f(x p ), (51) 
Jo Jo Jo m=1 J 



£(n,fc,x) := ^\x)^\l-x) 

Iff 1 \ . , , 2isin(nfc) / 1 



1 + Z2 ) cos ( nfc ) + ( 1 - ) cos[nfc(2a; - 1)] } . (52) 



In particular, 



rl 

Fi(n,k)= I dx£(n,k,x)f(x). (53) 



In the remainder of this section we explore the application of perturbation theory for treating the following equation 
whose real solutions yield the spectral singularities. 

oo 

ri-(fc) = 3^ = 0. ( 54 ) 

Suppose that we have a spectral singularity for the unperturbed potential, i.e., when 32 = 0. Let (rto,fco) be the 
value of (n,k) at which this spectral singularity is realized. If we turn on the perturbation, i.e., set ej2 7^ 0, this 
spectral singularity occurs for a new value (n*, fc*) of (n, k). Our aim is to express fc* and n* as power series in the 
perturbation parameter e, 

00 00 

m— m— 

and determine the coefficients k m and n m , that respectively take real and complex values. To do this we expand 
Fi(n, k) in a power series about (tto, &o)? 

n 

Ft(n, k)=J2 F ^( n ~ n o) p (k - k )\ (56) 

p,q=0 



7 




FIG. 1: (Color online) Schematic view of an infinite planar slab of gain material of thickness L that is aligned in the x-y plane. 



substitute k = fc* and n = ru in (|54|) . and use (|55[) and ([56| to express the resulting equation in the form 

oo 

5>e»' = 0, (58) 

3=1 

where cj are complex coefficients depending on k m , n m , and 32- For example, 

ci = Foiom + F 00 i fci +F10032, (59) 
C2 = ^010 «2 + F 00 i k 2 + -F020 txf + i^ou tiifci + -F002 k{ 

+ (^110 111+^101^)32+^200 32- (60) 

Finally, we enforce (|58j) by demanding that Cj = for all j > 1. Because Cj involves fc m and n m with to < j, in this 
way we obtain an infinite set of algebraic equations for k m and n m that we can solve iteratively. 
For example, in view of (|59")l and Ci = and C2 = give 

-Fbio n i + ^boi k\ = — -Fi.,0,0 3a s (61) 
Fqw n 2 + -fooi &2 = — (-P020 n i + -foil n i^'i + -F002 fc?) 

-(Fii ni- J Pioifci)32-i ? 20o3i (62) 

respectively. In general if n and k are independent complex and real variables, we can choose k m = and use the 
above method to compute n m . However, as we see in the next section, there are situations of physical interest where 
n depends on k and another real variable g. In this case, our method produces a perturbative calculation of k and g. 



V. OSSS OF AN INHOMOGENEOUS PLANAR SLAB GAIN MEDIUM 

A simple model supporting OSSs is an infinite planar slab gain medium [!|4|. Suppose we choose a coordinate 
system in which the slab is aligned parallel to the x-y plane and has a thickness L, as demonstrated in Figure[TJ Then 
the wave equation associated with this system admits a solution that propagates along the z-axis and corresponds to 
the following expression for the electric field ■ 



E(z,t) = E e~' luJt ^(z)e : 
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where Eq is a constant coefficient, uj is the angular frequency of the wave, e x is the unit vector pointing along the 
positive x-axis, ^ is a solution of the Schrodinger equation, 

-*"(*) + V(*)*(*) = ^*(«), (63) 
lj 2 \1 — n(u), z) 2 } 

— t , for z < L 2, 

V(z) := { c (64) 

for \z\ > L/2, 

c is the speed of light in vacuum, and n(u>, z) is the complex refractive index of the medium. As shown in Ref. 0, by 
a simple change of variables z and u>, namely 

z 1 Luj 

z^x :=- + -, uj^k:= , (65) 

LI c 

we can map (|63[) and (|64[) to |T]) and ([5]), respectively. This means that we can use our general results to compute 
the effects of inhomogeneity of the gain medium (z-dependence of the refractive index) provided that we express the 
perturbation parameter e, the coupling constants 31 and 32, and the function / entering (|44j) in terms of the physical 
parameters of the system. This requires the knowledge of the dispersion relation that determines the dependence of 
the complex refractive index on u and z. 

If we consider a gain medium that is obtained by doping a host medium of refraction index no and modeled as 
a two- level atomic system with lower and upper level population densities Ni and N u , resonance frequency u)q, and 
damping coefficient 7, then 

n*(u>,z) = n%- , 2 "ff., , (66) 
ur — 1 + 17 oj 

where Cj := lo/ujq, 7 := y/uio, Up ■— (iVj — iV u )e 2 /(m e eo), e is electron's charge, and m e is its mass. Furthermore, we 
have [I[ 

Mz) 2 = 2^{z)Jnl + n{z)^ 2jn n(z) = K (z) := (67) 

v Z7T 47T 

where Ao := 2hc/ujq is the resonance wavelength, g{z) is the effective gain coefficient (gain coefficient minus loss 
coefficient) at the resonance frequency, and we have made use of the fact that for all known gain media, k(z) <!C uq. 
Substituting (|67"1) in (|66)) gives 

n (u>, z) = n Q + — — — — — , (68 

2ir(ui z — 1 + 470;) 

Next, we assume that the pumping intensity decays exponentially inside the sample. This implies that if we pump 
it from the left-hand side, we have [7j, 

g(z) = (g + a) e -"(f+*> - a for \z\ < f, (69) 

where go := g( — %-), a is the attenuation coefficient at the resonance frequency that coincides with the largest allowed 
value of go, and v is the decay constant that specifies the exponential decay of the intensity of the pumping beam 
inside the slab. If we pump the sample from both sides (double pumping), we instead find [7j 



go + a 
cosh(f) 



cosh(^)-a for \z\ < 4. (70) 



Putting all these information together and making the change of variable (|55|) . we can reduce the problem of finding 
the OSSs of the above system to locating the spectral singularities of the potential ((6]) with vq and V\ given by (|44]) 
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and the following choices for 31, 32, e, and f(x), Q. 

3i 



( 2jtLu \' 
V A J 



1 - lln 



2ir 



32 :=2^ 2 f^ + l)3, 



a 



A 1 L 2 a A fnov for single pumping, 



A 1 L 2 a^fnoi /2 for double pumping, 



where we have introduced: 



cosh[^(x — i)] — cosh(^) 
v 2 cosh(|) 



for single pumping, 
for double pumping, 



(71) 
(72) 

(73) 



3 := 



1 



1 -Co 2 



We also note that according to (j47|) and (|7Tj) . 



n = \jnl 



7n A .go3 



7»^oAo3oA 2 



2- 



2tt(A 2 -r/A A- A 2 )' 



where A stands for the wavelength of the wave, i.e., 



A := 



2nc 2itL Aq 



(74) 



(75) 



(76) 



ui k Cj 

Because the physical parameters of practical interest are the wavelength A and the gain coefficient go , we fix all 
the other physical quantities and study the effect of inhomogeneity of the medium on the A and go values associated 
with spectral singularities. These we respectively denote by A* and g* and expand in power series in the perturbation 
parameter e, 



A* = A 



(0) 



5> 



9* = .9(0) + 9meT 



(77) 



Here A(o) and <?(o) are respectively the values of A and go associated with the spectral singularity of the unperturbed 
potential that appears for n = no and k = fco- As shown in Ref. Q, these are labeled by a mode number m and can 
be computed by substituting (|68|) and Co = A /A in (|49| and finding the real values of A and go that satisfy (|49|). In 
view of (|55|) and (|76|). it is easy to see that 



k 



2ttL 



A 



fci 



(0) 



2ttLAi 

A (o) 



fc Ai 



A 



(0) 



(78) 



Next, we view n as a function of A and go, and identify it with its Taylor series about (A(o), ff(o))- hi light of (|77|) . this 
gives for the value of n at (A*,<7*) the second equation in (|55t . i.e., n* = ^"^nme" 1 , and allows for identifying the 
coefficients n m in terms of those of the power series (|77|) for A* and g±. In particular, no is the value of n at (A(o) , <?(o)), 
i.e., 



3(i) 



3i 



^— \o) '9o—9(o) 



/2itLY 

vw 



"o-l- 



7?T,o A A 2 0) g ( o) 



27r ( A (o) ~ «7 A oA( ) - Aq) 



and ni is given by 



ni := ni, Ai + n ,iffi, 



dn 



«o,i 



dn 

dgo 



^— ^(0) > 9o—9(o) 



7»qA 2 A ( o)(2A + i7A(o))ff( ) 
47rn (A 2 - i 7 A A (0 ) - A§) 2 



-7?ioA A 2 0) 



47rn (A 2 0) - i7^oA (0 ) - Ag) 



(79) 
(80) 

(81) 
(82) 

(83) 
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Now, we recall that the spectral singularities are obtained from the equation 

oo 

r!- = X> 3 ^ = o. 



(84) 



In order to use perturbation theory to solve this equation for A and go, we need to express Fg and 32 & s functions of 
A, go, and e, substitute their Taylor series expansion about (A(o),3(o)) m 6 ) an d then set the coefficients 

of the resulting powers series to zero. As we will see below, for most realistic situations the first-order perturbation 
theory gives highly reliable results. Therefore, we outline the details of the calculation of the first-order corrections 
to A(o) and gro)> namely Ai and g\. 

Determination of Ai and g\ requires the following expansions of Fq, F\, and 32. 



Fo = F ( o) + (XX! + Ygi)e + 0(e 2 ), 



Ft 



F [1)+ 0{e l ), 32=3( 2) +0(e 1 ), 



where 0(e ) stands for the terms of order e and higher, and 



ip ip / 1 \ ip v ip 27tLF oi 
^(0) : = ro\?-o,Ko) = torn, A := nio-roio T2 



x (o) 



Y 

3(2) 



noi-Fbio 
32 



F(i) '■= Fi(n 0l ko) — F100, 
_ 27r -^o(.9(o) +«) 

^=^(0) I 90=9(0) 



"( A (0) 



«7AoA (0) - XI) ' 



(85) 

(86) 
(87) 
(88) 



In the derivation of these formulas we have made use ([55]). (|57)l. ([7T]) , ([?4"]) , ([75]). and (|8"T|) . 

Next, we observe that because for (n, k) = (no, ko) we have a spectral singularity, F( ) = 0. In view of this relation, 
~~ }, and (1571) . we can calculate 



-F010 — 



3(i) + i2k 
n 



1 

n 



3(i) 



X 



(0) 



^001 = - 1) = 7^ = 

K 



A 



(o)3(i) 
2nL 



These together with ([86]) and (jSSJ give 

mo 



no 



3(i) 



47Tli 



A 



(0) 



3(i) 



X 



Y 



(0) 



npi 
no 



3(i) 



47riL 
A(o) 



Furthermore, in light of Fro) = and (|84|) - ([88]), we can write the equation determining Ai and g\ as 

XX\ +Ygi = -Fioo3(2)- 

This is a complex linear equation involving two real unknowns. Therefore, we can easily solve it to obtain: 

Im(F 100 3(2)^*) . Im(F 100 3(2)^*) 



A 1 = - 



Im(XY*) 



jm(XY*) 



(89) 



(90) 



(91) 



(92) 



where "Im(-)" stands for the imaginary part of its argument. 

Next, we examine -Fioo that encodes all the information about the inhomogencity of the gain medium. According 
to JED, 



■fioo = iq(n ,/c ) 



1 



cfe£(n , k ,x)f(x), 



(93) 



no k a Jo 

where £ is given by (|52[) . With the help of ([49]) we have been able to simplify the expression for £(no, ko, x) and find 

1 



£(n ,fco,x) 



1 



[n k (x - |)] 



(94) 



Substituting this equation and ([73]) in ([93]) . evaluating the resulting integrals, and using ([49]) to simplify the outcome 
we obtain 



-fioo — 



(1 - (1 - e-")[4i + fc (l + njj)] - k v + 2iv 3 - 4k 3 n 3 (n 2 - l){v + e~ v - 1) 



2fcg n 4 v 2 {±k 2 n 2 + v 2 ) 



(95) 
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FIG. 2: (Color online) Location of the optical spectral singularities considered in Table [T] for < is < 0.5. The left- and 
right-hand figures respectively show the behavior of OSS for the singly and doubly pumped samples. The displayed dots that 
mark the boundaries of each of the curves of the OSSs correspond to v = and v = 0.5. As we increase v the location of OSSs 
in the A-po plane moves upwards on the given curves. 



for a singly-pumped sample, and 

P _ (n -l)(n + l) 2 [-2 + n^-2)-nH 

^ 100_ 2fc 2 n^(4fc 2 n 2 + ' W 

for a doubly-pumped sample. 

Now, we can use (|72|) . ([80]) . (|90l) . (|92|) . (|95|) . and (|96|) to calculate the first-order corrections to the wavelength and 
gain coefficient of the OSSs, i.e., Aie and g\e. 

As a concrete example, consider a semiconductor gain medium with the following specifications that is also studied 
in Refs. 



n 



= 3.4, L = 300 Aim, A = 1500 ran, 7 = 0.02, a = 200 cm^ 1 , < v < 0.5. (97) 



For this sample e 0.8 < 0.4 for single pumping and e ~ 0.8 is 2 < 0.2 for double pumping. In particular for 
physically realistic situations where v ^ 0.1, we have e 2 ^ 6.4 x 10^ 3 and e 2 ^ 4.0 x 10 -4 for single and double 
pumping, respectively. This shows that the first-order perturbation theory produces a highly reliable description of 
the OSSs for this system. 

In the remainder of this section, we report the results of the first-order perturbative calculation of g+ and A* for 
the first five OSSs that appear as we increase the intensity of the pumping beam starting from zero. These have a 
wavelength that is closest to the resonance wavelength of the sample, Ao = 1500 nm. We label them using the mode 
number m that is introduced in Q and takes values between 1358 and 1362, with 1360 corresponding to the OSS 
with closest wavelength to Ao- Table U gives numerical values of A* and g± for these OSSs and five different values 
of the decay constant v. As we increase is, A* remains essentially unchanged while g* increases. This is particularly 
pronounced for a singly-pumped sample. It confirms the semiclassical results obtained in 0- 

Figure [5] shows the curves traced by (A*,<7*) as we change v for each of the OSSs considered in Table [J] The fact 
that these curved are essentially vertical line segments shows that A* docs not depend on is, while the opposite is 
true for g±, especially for the singly-pumped sample. The displayed dots that mark the lower and upper boundaries 
of each line segment respectively correspond to is = and is = 0.5. As we increase is, the location of each OSS in the 
A-(?o plane moves upwards along the corresponding line segment. 
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Single Pumping 


Double Pumping 


m 


V 


A (mn) 


go (cm -1 ) 


A (nm) 


go (cm -1 ) 




0.0 


1497.561770810 


41.53101 


1497.561770810 


41.53101 


1362 


0.1 


1497.561770785 


43.45261 


1497.561770784 


41.56407 


0.2 


1497.561770716 


45.25128 


1497.561770707 


41.66286 




0.3 


1497.561770609 


46.93581 


1497.561770579 


41.82620 




0.5 


1497.561770304 


49.99447 


1497.561770180 


42.33818 




0.0 


1498.389018373 


40.91032 


1498.389018373 


40.91032 


1361 


0.1 


1498.389018341 


42.83283 


1498.389018339 


40.94324 


0.2 


1498.389018251 


44.63206 


1498.389018239 


41.04159 




0.3 


1498.389018115 


46.31686 


1498.389018073 


41.20423 




0.5 


1498.389017715 


49.37529 


1498.389017554 


41.71399 




0.0 


1499.999983312 


40.40905 


1499.999983312 


40.40905 


1360 


0.1 


1499.999983275 


42.33379 


1499.999983220 


40.44217 


0.2 


1499.999983205 


44.13541 


1499.999983115 


40.54115 




0.3 


1499.999983098 


45.82273 


1499.999983003 


40.70480 




0.5 


1499.999982791 


48.88649 


1499.999982512 


41.21777 




0.0 


1501.475689102 


40.79650 


1501.475689102 


40.79650 


1359 


0.1 


1501.475689077 


42.72315 


1501.475689075 


40.82968 


0.2 


1501.475689007 


44.52660 


1501.475688997 


40.92881 




0.3 


1501.475688899 


46.21566 


1501.475688689 


41.09272 




0.5 


1501.475688590 


49.28266 


1501.475688464 


41.60649 




0.0 


1502.670951310 


41.63220 


1502.670951310 


41.63220 


1358 


0.1 


1502.670951286 


43.56043 


1502.670951282 


41.65321 


0.2 


1502.670951220 


45.36542 


1502.670951211 


41.76466 




0.3 


1502.670951118 


47.05600 


1502.670951089 


41.92890 




0.5 


1502.670950826 


50.12593 


1502.670950707 


42.44373 



TABLE I: The wavelength A* and gain coefficient for the five OSSs that are generated by pumping the semiconductor gain 
medium (|97|l . m is the mode number labeling these OSSs Q. v is the decay constant for the intensity of the pumping beam(s) 
inside the sample. 



VI. CONCLUDING REMARKS 

In this article we have outlined a general method for carrying out a perturbative calculation of the transfer matrix 
for a general complex potential that vanishes outside a closed interval. This allows for a systematic characterization 
of the spectral singularities of this class of potentials. This turns out to be particularly suitable for the study of 
optical spectral singularities and their time-reversed analogs that respectively correspond to lasing at threshold gain 
and antilasing. 

As an application of our general results, we examined the problem of computing the transfer matrix for a potential 
consisting of finitely many Dirac delta-functions that are centered at arbitrary points and have arbitrary complex 
coupling constants. We showed that perturbation theory gives an exact expression for the transfer matrix for this 
system. 

Next, we considered an arbitrary complex perturbation of a constant (complex) barrier potential and applied our 
method to compute the effect of the perturbation on the spectral singularities. A physical realization of this model 
is in the description of threshold lasing associated with an infinite planar slab of gain material. For this system 
the gain coefficient that is proportional to the intensity of the pumping beam decays exponentially as the beam 
penetrates the medium. This in turn makes the gain medium inhomogencous. Our method allows for an essentially 
analytic calculation of the effect of this inhomogeneity on the location of optical spectral singularities. Our results 
confirm those obtained using the method of Ref. [7| that is based on the semiclassical approximation. Compared with 
this method ours has the advantage of being applicable in every spectral range. In particular, we can use it in the 
spectral ranges comparable with the length scale of the system. This is, for example, the case in the recent study 
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of unidirectional invisibility jl8j . Our method allows for a thorough analysis of this and much more general optical 
systems displaying threshold lasing, antilasing, and unidirectional invisibility. 
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